Near-zero-energy end states in topologically trivial spin-orbit coupled superconducting nanowires 
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A one-dimensional spin-orbit coupled nanowire with proximity-induced pairing from a nearby s-wave super- 
conductor may be in a topological nontrivial state, in which it has a zero energy Majorana bound state at each 
end. We find that the topological trivial phase may have fermionic end states with an exponentially small energy, 
if the confinement potential at the wire's ends is smooth. The possible existence of such near-zero energy levels 
implies that the mere observation of a zero-bias peak in the tunneling conductance is not an exclusive signature 
of a topological superconducting phase even in the ideal clean single channel limit. 



In one dimension, topological superconducting wires are 
predicted to support a localized Majorana bound state at each 
end These Majorana states are particle-hole symmetric 
and have exactly zero excitation energy. Within the associated 
degenerate subspace, braiding and exchange operations can be 
shown to be non-Abelian |2-4|, making them potentially use- 
ful in a topological quantum computation schemes lISHH). It 
has been recently noted, that in the right parameter regime, 
spin-orbit coupled semiconductor nanowires with proximity 
induced superconductivity, should exhibit the required topo- 
logical superconductivity for Majorana pair formation 191 [TOl . 

In light of these proposals, the experimental observations of 
zero-bias peaks in normal-metal superconductor tunnel junc- 
tions, which are unaffected by small variations of the mag- 
netic field or gate voltages, may indicate the presence of topo- 
logical superconductivity ifTTl [T2l . While these observations 
are a necessary indicator of the predicted mid-gap Majorana 
states ifTJl [141 . it is crucial that alternative mechanisms for 
the zero-bias conductance be ruled out in order for them to 
be decisive. One example of such an alternative mechanism 
applies to quasi-one dimensional wires with multiple conduct- 
ing channels 1 15 -19|, for which low-energy fermionic bound 
states are predicted to appear in the topological as well as in 
the non-topological phase if the Zeeman energy exceeds the 
splitting between transverse subbands |20-24|. Other alterna- 
tive mechanisms involve disorder |25 , 26 1, possibly in combi- 
nation with a gapless region at the wire's end |27|. These latter 
findings suggest that clean single-channel wires offer a favor- 
able setting to discern the presence of Majorana end states. 
Indeed, experiments are progressively approaching this ideal 
scenario fTTlT^J. 

In this letter we show that the original proposals iP.TOl for 
topological superconductivity in clean one-dimensional semi- 
conductor wires also allow for near-zero-energy end states 
deep in the topologically trivial phase, provided the potential 
that confines the electrons at the wire's end is smooth. The low 
energy is a systematic property of these states, that persists as 
long as the confining potential and the induced superconduc- 
tivity are smooth functions of position. The existence of such 
low-energy Andreev states leads to a low-energy peak in the 
tunneling conductance in the topologically trivial phase. Since 
gate-induced confinement potentials are typically smooth, the 
mechanism we describe here may be relevant for the recent 



experiments lfm[T2l . Our analysis is consistent with and ex- 
plains the observation of zero-bias conductance peaks in re- 
cent numerical simulations of clean semiconductor wires by 
Pradaefa/. 128)1 . 

Following the original theoretical proposals |9 10|, we 
consider a one-dimensional semiconductor with Rashba spin- 
orbit coupling of strength a, subject to a magnetic field with 
Zeeman energy B > and proximity-coupled to a standard s- 
wave spin-singlet superconductor. Such a system is described 
by the four-component Bogoliubov-de Gennes Hamiltonian 

H = + V{x) -II- Ba^ + apay^ + AayT^, (1) 

where (T^^y^z and r^.y^z are Pauli matrices acting on the spin 
and particle-hole degrees of freedom, respectively. Further, 
m is the effective electron mass, /i — pp/2m the chemical 
potential, A the proximity induced superconducting gap in the 
absence of the magnetic field, and V{x) is the potential that 
describes the confinement of electrons near the wire's end. 

As shown in Refs. |9| and [TOl the Hamiltonian ([TJ is in a 
topological phase with Majorana fermions at its ends if i? > 
Be — \/ fi^ + A^. Here we consider the topologically trivial 
regime with weak induced superconductivity, A ^ B ^ /i. 
The condition S 3> A rules out spin singlet s-wave pairing, 
so that the induced superconductivity must be of p-wave type. 
However, unlike in the topological regime, where the model 
([T]i effectively admits p-wave superconductivity for one spin 
channel only, in the non-topological regime _B ^ /i both spin 
channels acquire superconducting correlations. If B ^ A the 
two spin channels exist as effectively independent p-wave su- 
perconductors in the wire's bulk, but they are coupled at the 
wire's ends, which gaps out the pair of Majorana-like excita- 
tions that would have existed at the wire's end for uncoupled 
channels. As we show below, this coupling is strong if the 
wire's end is abrupt, but weak if the confinement is smooth, 
which explains the appearance of an Andreev bound state at 
an energy far below the bulk excitation gap. The crucial dif- 
ference between an abrupt ending and a smooth confinement 
is that the spin-orbit energy Ego — ap remains finite up to the 
turning point for a hard-wall confinement, whereas Eso goes 
to zero continously for a smooth confinement. 

In order to arrive at an approximate analytical solution of 
this problem, we assume that the energies B, Eso, and A, are 
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much smaller than the kinetic energy /i — V{x). This separa- 
tion of energy scales breaks down near the turning point at the 
wire's end, where the velocity 



v{x) = y/2[fi-V{x)]/m 



(2) 



goes to zero. We circumvent this difficulty by solving a mod- 
ified version of the problem, in which the wire has a hard- 
wall confinement with V = inside the wire, and a position- 
dependent spin-orbit strength a with 



£soix) ~ mav{x) — a{x)pi 



(3) 



to account for the position-dependence of Ego- The two de- 
scriptions are essentially equivalent if the coordinate x in the 
hard-wall model is related to the original coordinate x as 



dx 



v(x) ' 



(4) 



where up ~ ■\/2/i/m is the Fermi velocity in the hard-wall 
model. Taken together, the relations ([3]l and Q ensure that 
the electrons "see" the same Zeeman energy B and spin-orbit 
energy Ego as a function of time when they reverse their direc- 
tion at the wire's end. 

The inequalities B, Eso, A ^ /i allow us to linearize the 
kinetic energy, writing 



(5) 



where the functions ip± are slow functions of position on the 
scale h/pp. The function -0+ describes right-moving elec- 
trons and left-moving holes, while describes left-moving 
electrons and right-moving holes. They are subject to the four- 
component Bogoliubov-de Gennes Hamiltonian 

H± = TifiVFTzdi - Ba^Tz ± a{x)p-payT:, + AayT^,(6) 

and the boundary condition ^+{0} = e^*'''i/'-(0), rj being a 
phase shift characteristic of the detailed boundary conditions 
at the wire's end at i = 0. 

The normal part of the Hamiltonian (|6| can be diagonalized 
by a rotation in spin space. Defining the angle 6{x) and the 
wavenumber fcni(i) > as 



B = hvFk^{x) cos 0{x), 
£soiS:) — hvpk^{x) sm9{x), 



(7) 



a basis change maps the Bogoliubov-de Gennes Hamiltonian 
(jejl to ffo,± + Hi_±, with 

Hq,± = hvp{^idx - km<7z)Tz ± A<JzTx sin6'(i), (8) 



Hl^± = AOyTx cos,d{x) 



hvp 89 



(9) 



The superconducting pairing in Hq ± is of p-wave type and 
pairs electrons of equal spin (in the rotated frame) with p- 
wave gap 



Ap = A sin e 



^soA 



(10) 



Analytical 
Tight binding - 




FIG. 1. (Color online) Andreev end-state energy e, normalized to the 
bulk excitation gap Ap, as a function of adiabaticity parameter a. 
The red squares are obtained from a numerical calculation for which 
the wire is terminated by a smooth potential V{x) of the fonn 1 17 i 
with a = 5/^. Other parameters used in this numerical calculation 
are: £so = 0.1/^, A = 0.05/i, B — 0.275^. The solid blue curve 
corresponds to Eq. ( |13[ l of the main text. 



whereas the superconducting pairing in Hi ± is of s-wave 
type and connects electrons of opposite spin. 

In the limit B ^ A and for a smooth confining potential, 
Hi,± can be treated in perturbation theory. The unperturbed 
Hamiltonian -ffo.± admits two zero-energy end states of Ma- 
jorana type. 
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where the superconducting coherence length ^ is defined as 
hvp/£^{x) — A| sin6'(x)| and $7 is a normalization constant. 
Calculating the matrix element ofHi ± between these states, 
we find that the wire's end harbors a single Andreev end state 
with energy 



X cos 



dx 



do 

2Acos6'(i) + hvp — 
ox 



dx'ktn{x') 



(12) 



Returning to the parameters of the original model ([T]|, the en- 
ergy e of the Andreev end state reads 



IB 



°° , 2AJWTI^Jxf ~ haidV I dx) 

nT 

v{x)[B^ + eU^r] 



X cos 



2 / dx'k^{x') 

Xg 



^-2f:^dx'i/ax') 



(13) 
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where Hv{x)k^{x) = ma(,{x)A/H = B"^ + e^^ and 



= 4 / dx- 



v{x) 



(14) 



For the simple example that V{x) has a linear dependence on 
x near the wire's end, y(a;) = /x—y'a;, with the condition that 
V ^ BA/ha — which ensures that that Ego <C B through- 
out the entire integration range — , a closed-form expression 
can be obtained and one finds 



(15) 



This result should be compared with the energy of the An- 
dreev end state for a hard wall, which is 



(16) 



if Eso, A <^ B. The energy ( |T5) l is essentially zero — even 
if compared with the p-wave gap Ap — for a range of mag- 
netic field far below the critical field Be at which the transition 
to the topological phase takes place. (Incidentally, even with 
hard-wall boundary conditions, the Andreev end-state energy 
£ may be small in comparison to Ap if B is sufficiently large 
in comparison to Sso and A.) Experimentally, the difference 



between the finite excitation energy e of Eq. (15 1 and the strict 
zero-energy of the Majorana bound states may be difficult to 
resolve. 

Both the splitting from the singlet pairing [first term in Eq. 
(j9]l] and the splitting from the non-adiabaticity of the confin- 
ing potential [second term in Eq. (j9]l] vanish exponentially in 
the limit of a smooth, adiabatic confinement, provided the 
Zeeman energy B sufficiently far exceeds A. This is the 
main result of this letter. We stress that in the weak pair- 
ing limit, the splitting e decreases with increased coherence 
length, and therefore cannot be simply understood as the re- 
sult of the (small) separation between the turning points for 
electrons with opposite spin. (Indeed, no such separation is 
present in the effective hard-wall model used for our calcula- 
tion!) 

We have compared the theoretical predictions to numeri- 
cal tight-binding simulations of a discretized version of the 
Hamiltonian (jTJ. For this purpose we choose the confining 
potential 



Vix) = ( -"^^''"^^ - > 0' 
^ ' [oo X <0, 



(17) 



where the parameter a controls the degree of adiabaticity. Fig- 
ures [T] and |2] show representative results of the numerical cal- 
culation, together with the analytical prediction of Eq. 
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As discussed above, the small energy e of the Andreev end- 
states results from the ineffectiveness of a smooth potential to 
couple the two Majorana modes for the two spin channels. 
This near-degeneracy will be lifted in the presence of per- 
turbations with an abrupt spatial dependence that couple the 
different spin-orbit bands. Examples of such perturbations 
are scattering from point-like impurities (which couple left- 
moving and right-moving particles), or a the abrupt vanish- 
ing of the pairing potential, which happens, e.g., if not all of 
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FIG. 2. (Color online) Andreev end-state energy e as a function of 
the Zeeman energy B (main panel) and induced superconducting gap 
parameter A (inset). The parameters of the numerical calculation 
are: Ebo = O.l/i, A = 0.04/^ (main figure), B = 0.5^ (inset). 
The confining potential has the form (jTTJ with a = 5fi and a = 
1.273(/i/pp). 



the semiconducting wire is covered with the superconducting 
contact. The Andreev end-state energy e in the presence of a 
point impurity with potential US{x — Xi)Tz is (to first order in 

U) 



iUame 



n^B^ + eU^iy 



sm 



4?7 



dx 



mv(x) 



(i8) 



For the example of a slowly varying potential V{x) = fi—xV 
with a linear dependence on position, this gives 



aWA 



- [Xi-Xo) 



4??+ 1^ {x.-xofl^ 



(19) 



In the case that the order parameter vanishes abruptly. 



3J1, the end-state en- 
T3]l. For the special 



A{x) = AQ{x — xn) (see inset of Fig 
ergy can be obtained directly from Eq. 
case that the entire potential modulation occurs in the normal 
region, the discontinuity in A contributes to the end-state en- 
ergy by the amount 



2hB 1 

Q, mavp 



-Re 



2i dxk^Jx) 

1 + ik„,£, 



where 



xn 



4 / dx—- 

'xo 



(20) 



and fcm, ^, and vp are the asymptotic values for x > x^. Fig- 
ure [3] compares numerical simulations of the model ([T]) with 
and without an abrupt change in the superconducting order 
parameter A. 

Up to this point, our discussion has focused on one- 
dimensional semiconductor wires with a single transverse 
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FIG. 3. (Color online) Andreev end-state energy e versus the in- 
duced superconducting gap parameter A with (squares) and without 
(circles) an abrupt termination of the order parameter A at position 
a;N. The inset shows the functional form of the superconducting or- 
der A{x) = AO{x — xn). In both cases the wire is terminated by a 
smooth potential V{x) of the form 1 17 1 with parameters with a = Sfi 
and a ~ 3.183/i/pf. Other parameters used in this numerical calcu- 
lation are e^o = O.l/i and B = 0.275^. 



channel. Our arguments continue to be valid for multichan- 
nel wires. In this case, for B > A each transverse channel 
is in a separate effectively spinless p-wave superconducting 
state. Hard-wall boundary conditions at the wire's end cou- 
ple the channels, which gaps out the end states, up to the 
possible exception of a single Majorana end state if the to- 
tal number of channels (counting spin) is odd. In general, 
a coupling exists between spin degenerate channels with the 
same transverse mode, as well as between channels with dif- 
ferent transverse modes, although the "off-diagonal" coupling 
is small if the wire width is much smaller than the supercon- 
ducting coherence length ^ because of an approximate chi- 
ral symmetry ll20l|29i . By the mechanism discussed above, a 
smooth confinement strongly reduces the diagonal and the off- 
diagonal couplings between channels, giving rise to int {N/2) 
low-energy fermionic states at the wire's end. This is illus- 
trated in Fig. |4]for the case of a semiconducting wire with two 
spin-degenerate transverse channels. 

In conclusion, we have shown that a sufficiently smooth 



confinement potential at the wire's ends leads to the exis- 
tence of low-energy Andreev end states even in the topo- 
logically trivial phase of one dimensional proximity coupled 
nanowires. These results could be relevant for recent experi- 
ments lfTn[r2 l . in which the confinement at the nanowire ends 
is gate-induced. The presence of such low energy Andreev 
states would lead to near zero conductance peaks deep in the 
topologically trivial parameter regime. This low energy peak 
hinders the experimental verification of a topological super- 
conducting phase via tunneling density of states, even in the 
ideal single channel case. 

We gratefully acknowledge discussions with J. Danon, F. 
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by the Alexander von Humboldt Foundation in the framework 
of the Alexander von Humboldt Professorship, endowed by 
the Federal Ministry of Education and Research. 
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FIG. 4. (Color online) Andreev end-state energies e for a two- 
dimensional wire with two transverse channels, normalized by the 
bulk excitation gap, as a function of the adiabaticity parameter a. 
The Hamiltonian is given by the two-dimensional extension of Eq. 
l[T](, which has the spin-orbit coupling term apxO-yTz — apyOx. Pa- 
rameters in the numerical calculation are B = 0.1667/i, and wire 
width W = 1.225/i/pf = 0.054^. The spin-orbit energies of the 
two transverse bands (calculated for B — 0) are Ebo.i ~ 0.074/^, 
£bo,2 = 0.042^1. 
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